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Abstract. We develop an analytic model for galaxy intrinsic alignments (IA) based on the 
theory of tidal alignment. We calculate all relevant nonlinear corrections at one-loop order, 
including effects from nonlinear density evolution, galaxy biasing, and source density weight¬ 
ing. Contributions from density weighting are found to be particularly important and lead 
to bias dependence of the IA amplitude, even on large scales. This effect may be responsible 
for much of the luminosity dependence in IA observations. The increase in IA amplitude 
for more highly biased galaxies reflects their locations in regions with large tidal fields. We 
also consider the impact of smoothing the tidal field on halo scales. We compare the perfor¬ 
mance of this consistent nonlinear model in describing the observed alignment of luminous 
red galaxies with the linear model as well as the frequently used “nonlinear alignment model,” 
finding a significant improvement on small and intermediate scales. We also show that the 
cross-correlation between density and IA (the “GI” term) can be effectively separated into 
source alignment and source clustering, and we accurately model the observed alignment 
down to the one-halo regime using the tidal field from the fully nonlinear halo-matter cross 
correlation. Inside the one-halo regime, the average alignment of galaxies with density trac¬ 
ers no longer follows the tidal alignment prediction, likely reflecting nonlinear processes that 
must be considered when modeling IA on these scales. Finally, we discuss tidal alignment in 
the context of cosmic shear measurements. 
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1 Introduction 

Coherent, large-scale correlations of the intrinsic shapes and orientations of galaxies are 
a potentially significant source of systematic error in weak gravitational lensing studies. 
Contamination from these galaxy “intrinsic alignments” (IA) can bias or degrade lensing 
science results. Since early work establishing its potential effects [1—4], IA has been examined 
through observations [5-8], analytic modeling [3, 9, 10], and simulations [11, 12] - see [13] 
for a recent review. As our understanding has improved, IA has also emerged as a potential 
probe of large-scale structure as well as halo and galaxy formation and evolution [14, 15]. 
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Despite recent progress, further understanding of IA is needed. The increased statistical 
precision of current and next-generation lensing surveys such as DES 1 , KiDS 2 , HSC 3 , Euclid 4 , 
LSST 5 , and WFIRST 6 will require unprecedented control of systematic uncertainties. It is 
likely that a combination of approaches will be employed for effective IA mitigation. We 
focus in this work on analytic modeling of galaxy alignments, usually relating IA to the 
gravitational tidal held, which is determined by the surrounding large-scale structure [3, 
9]. These analytic approaches are naturally perturbative and are analogous to galaxy bias, 
which relates the galaxy and dark matter density fields. On sufficiently large scales, where 
perturbative treatment of the tidal held converges, these models should produce a simple 
and accurate description of IA. On smaller scales, however, a perturbative description will 
necessarily break down, and a qualitatively nonlinear treatment is required. To model IA 
on these scales, halo model prescriptions have been developed to specify the distribution of 
shapes and orientations for central and satellite galaxies (e.g. [16, 17]). These halo models 
can also be combined with effective tidal held models on larger scales. 

Given the current challenges and uncertainties in modeling IA, especially at small 
scales and for the faint sources that dominate the lensing signal, it is advisable to also 
consider model-independent methods to mitigate IA contamination. However, more general 
approaches to IA mitigation tend to come at a signihcant cost of statical signal (e.g. [18]). 
For IA marginalization, it is important to reduce the number of “nuisance” parameters as 
well as accounting for the covariance of these parameters for multiple cosmological probes. 
Well-understood and accurate analytic modeling allows for efficient parameterization of IA 
for use in large joint analyses. 

First proposed in [3], the tidal alignment model posits that halo and galaxy shapes (or 
the correlated components) are aligned with the local tidal held. This type of alignment 
is most likely to arise in elliptical galaxies for which angular momentum does not play a 
signihcant role in determining shape and orientation. Indeed, tidal alignment has been shown 
to agree with shape correlations on large scales for luminous red galaxies (LRGs), which are 
highly biased, elliptical, and pressure-supported [6, 8, 10]. For spiral galaxies, the acquisition 
of angular momentum, for instance through “tidal torquing,” generally leads to quadratic 
dependence on the tidal held (e.g., [9, 19]). While tidal alignment may accurately capture 
the relevant astrophysics for certain types of galaxies, it can also be motivated through 
symmetry arguments. Since the tidal alignment contribution is the lowest-order function 
of the gravitational potential with the necessary symmetry, it is expected to dominate IA 
correlations on sufficiently large scales. Due to its simplicity and its relative success at 
describing the IA signal of highly biased, red galaxies (the only sample with well-measured 
shape alignments), the tidal alignment model remains an important tool. 

Because of the linear relationship assumed between galaxy shape and the tidal held, tidal 
alignment has historically been referred to as the “linear alignment” model. However, it was 
quickly recognized that nonlinear effects are likely important on smaller scales. Attempts 
to include nonlinear dark matter clustering by using a nonlinear matter power spectrum 
yielded the “non-linear alignment” (NLA) model, which provides qualitative improvement in 


1 Dark Energy Survey, https://www.darkenergysurvey.org 

2 Kilo Degree Survey, http://kids.strw.leidenuniv.nl/ 

3 Hyper Suprime-Cam, http://www.subarutelescope.org/Projects/HSC/ 

4 http://sci.esa.int/science-e/www/area/index.cfm?fareaid=102 

’Large Synoptic Survey Telescope, http://www.lsst.org 

®Wide-Field Infrared Survey Telescope, http://wfirst.gsfc.nasa.gov/ 
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matching LRG alignment on intermediate scales [9, 20]. The NLA model has been the fiducial 
functional form for several weak lensing projections and analyses (e.g. [6, 21]). However, 
once effects beyond linear order are included, referring to the “(non)linear alignment model” 
becomes ambiguous. Instead, we use the term “tidal alignment” to refer to any model in 
which the relationship between intrinsic galaxy shapes and the tidal field is linear. Quadratic 
alignment (e.g. tidal torquing) models involve a quadratic dependence on the tidal field. 

Within the context of tidal alignment, three effects can produce nonlinearities in the 
intrinsic shape correlations: (1) nonlinear evolution of the dark matter density field, leading 
to nonlinear evolution in the tidal field; (2) observing the IA field at the positions of shape 
tracers, leading to a weighing by the tracer density; (3) a nonlinear bias relationship between 
the galaxy and dark matter density fields. The NLA model includes the nonlinear evolution 
of the dark matter density but but does not consider other nonlinear contributions. While 
the NLA approach may improve the model fit to data, it is not fully consistent and omits 
important astrophysical effects. Some work (e.g. [22]) advocates explicitly ignoring the non¬ 
linear evolution of the density field when calculating the intrinsic ellipticity field, assuming 
that alignment is set at early times when the density field is closer to linear. As discussed 
below, this physical picture may be relevant, and the redshift at which IA is set is an impor¬ 
tant element in any model. In this work we wish to treat all nonlinear contributions on equal 
footing, providing a more complete description of the tidal alignment scenario. With such 
a framework in place, we can more convincingly examine how well tidal alignment captures 
the physics of IA and on what scales a nonlinear relationship between IA and the tidal field 
must be considered. Having a more complete treatment of possible nonlinear effects will also 
allow the extension of the tidal alignment model to the types of lensing sources that are 
actually relevant for current and upcoming studies. We now have several measurements of 
IA for LRGs, but these are not typical lensing sources. Thus, even if the ad-hoc NLA model 
is satisfactory for LRGs, it may not be for the fainter sources that will dominate the lensing 
signal. 

We employ two distinct approaches to examine nonlinear tidal alignment. First, we 
use a standard perturbation theory (SPT) expansion of the density field to calculate the 
one-loop non-linear contributions to the tidal alignment model. Within this expansion, we 
examine the effect of smoothing the tidal field, an element of IA modeling which remains 
an important open question (e.g. [10, 15]). As a second approach, we show that source 
clustering can be effectively separated from tidal alignment itself, allowing shape correlations 
to be described as a product of the two effects. We then employ non-perturbative modeling 
to describe the tidal field around a dark matter halo, allowing a description of IA down 
to significantly smaller scales. In the one-halo regime, we find that recent IA measurements 
indicate relatively constant galaxy alignment as a function of separation, indicating processes 
beyond tidal alignment. 

This paper is organized as follows. Section 2 outlines the relevant formalism for measur¬ 
ing shape correlations. Section 3 develops the tidal alignment model, consistently treating all 
nonlinear contributions using SPT, and describes the impact of smoothing the tidal field. We 
compare the model predictions with IA measurements. Section 4 explores a non-perturbative 
description of tidal alignment, separating the observed intrinsic shape correlations into source 
alignment and clustering terms. We also discuss the impact of nonlinear tidal alignment con¬ 
tributions in the context of cosmic shear. A discussion of the main results is in Section 5. 
The details of many of the relevant calculations are presented in an appendix. We assume a 
flat ACDM cosmology with Q m = 0.279, Qj, = 0.046, erg = 0.808, and h = 0.701. 
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2 Formalism of intrinsic shape correlations 


Observed shapes can be expressed in terms of an ellipticity, eo and position angle, 0, measured 
with respect to a chosen direction. The shape can also be decomposed into components 
(e_|_,e x ). In terms of the position angle and minor-to-major axis ratio, b/a: 


e+ 

e x 


(l-(b/a) 2 \ 

cos(2</>) 


cos(2 4>) 

\1 + (b/a) 2 ) 

sin(20) 

— eo 

sin(2 4>) 


( 2 . 1 ) 


Note that this decomposition obeys the necessary spin-2 symmetry (invariance under rotation 
by nn). 

In the weak lensing regime, an ensemble of observed galaxy shapes provides an estimate 
of the lensing shear ( 7 ). The two quantities are related through the shear responsivity 1Z, 
which captures the average response of measured ellipticity to a small shear [23]: 

(7+>7x) = ^((e+,e x ). (2.2) 

The observed (projected) shape of a galaxy has contributions from both intrinsic shape (y 7 ) 
and gravitational shear (y G ). Since we are concerned with coherent alignment effects that 
will contribute to the ensemble shape average, in the remainder of this work, we will use 
shear and ellipticity interchangeably, denoting both with 7 . For small lensing shears, these 
contributions add: 7 ° bs = y G + j 1 . The intrinsic galaxy shapes will in general include both 
a random “shape noise” and a coherent intrinsic alignment. The lensing contribution to the 
observed shape is typically only ~ 1%. 

In configuration space, the galaxy density-weighted shape field (see section 3.1.2) is the 
product of the galaxy density and intrinsic shape fields: 7 (+, x ) = (l + 5 s )y(+ )X ), where we use 
the standard definition of the overdensity field: 5 = p/p — 1. In Fourier space, the natural 
decomposition of this shape field is into E (curl-free) and B (divergence-free) modes. E 
and B modes can be thought of as measuring the + and x components, with respect to the 
direction of k (e.g. [24]). 

7c( k ) = M k)7+(k) + /s(k)7 X (k), 

IB (k) = -/s(k)7+(k) + / E (k)7 x (k), (2.3) 

where we have defined the angular operators: 


Mk) 
Mk) 



‘Ikxky 


(2.4) 


for k 2 = k 2 + k 2 . We work in the plane-parallel approximation and define the z-direction 
to be along the line-of-sight. Because the decomposition into E and B modes is invariant 
under rotations on the plane of the sky, we can consider k y = 0 modes only, without loss of 
generality. In terms of these Fourier space fields and the galaxy density field <5 fl (k), we can 
define the following power spectra: P 9 e( k), Pee{ k), and Pbb{ k). By symmetry under parity, 
P g B( k) and Peb( k) must be zero. Similarly, the configuration space correlation functions 
£gx = C+x = 0. 
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From Equation 2.3, the relevant power spectra are: 


(27r) 3 i5(k + k.')P gE (k, 9 k ) = <«5 fl (k)fe(k')) = <<5 9 (k)|7 + (k')) fcl/= o, 

( 2ir) 3 5(k + k')P EE (k,d k ) = (^(kJl^k')) = (7 + (k)|7 + (k , )) fcj) =o, 

(2Tr) 3 5(k + 'k')PBB(k,e k ) = (7 jB (k)|7 B (k / )) = (7x (k)l7x (k'))fc y =o, (2.5) 


where the power spectra depend on the amplitude k as well as 9 k . the angle between k and 
the line-of-sight. In general, the two galaxy samples being correlated need not be same, 
particularly in the case of density-shape cross-correlations (e.g., [8]) 

Correlation functions in configuration space are defined for galaxy shape components 
measured with respect to the projected separation axis between pairs of objects [24], In 
terms of the power spectra defined in eq. 2.5: 

/ d 3 k 

j^^e^ r P gE {k,d k )cos(2(j) k ), (2.6) 

/ d 3 k 

^-^e* k ' r ( P E E(k , 8 k ) cos 2 {2(j) k ) + P B B(k, 9 k ) sin 2 (2 fa)) , 

/ d 3 k 

^^3 e 4k ' r ( P EE (k , d k ) sin 2 (2 <f> k ) + P B Ei(k, 0 k ) cos 2 (2 fa)) , 

where r p is projected (transverse) separation, and II is line-of-sight separation. Switching 
to cylindrical coordinates and performing the <p k part of the integral yields the standard 
relations in terms of Bessel functions J n (x): 


-2 


ig+{r pi n) = J dk z dn ^cos ( II) J 2 ( kt p )P gE (k), 


(2.7) 


£(++,xx)(rp,n) = ^ J dk z dnKCOs(k z U) 


x [(Jo(nr p ) ± Ji{nr p )) P EE (k,0 k ) + (J 0 (nr p ) J 4 (k7j,)) P B B(k, 0 k )\ . 


The wavevector has been decomposed into components along the line-of-sight [k z = kcos9 k ) 
and on the sky (k = ksin9 k ). The top line of eq. 2.7 uses P gE (k,9 k ) = sin 2 (9 k )P gE (k). 
The ^-dependence for P EE and P B b is more complicated, and this dependence is left in the 
definitions of the power spectra. 

Ilmax 

The projected correlation functions are defined as w(r p ) = f dn£(r p ,II), yielding: 

Ilmax 
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For separations r p <C n max , only k z = 0 modes contribute (the Limber approximation): 


OO 

-2n f 

Wg+(r p ) « J dKKj 2 (Kr p )Pg E {K), 


(2.9) 


7r / 

u)( ++iXX )(r p ) « ^ J dun [(Jo(Krp) ± J 4 (rcr p )) P E e(k) + {Jo(nr p ) J 4 («r p )) Pss(«)] • 


This result can be seen from taking the limit n max —>• oo in eq. 2.8. 

The E and B mode correlation functions can be written as the Fourier transform of the 
corresponding power spectrum: 


C ( EE , BB ) 


d 3 k 

(27t) 3 


e* k r P(EE,BB)(k). 


( 2 . 10 ) 


Equivalently, these correlations can be written as non-local functions of £ ++ and £ xx [10, 25]. 


3 Tidal alignment and nonlinear contributions 

The “tidal alignment” model (sometimes referred to as “linear alignment”) posits that in¬ 
trinsic galaxy ellipticity has a linear dependence on the tidal held [9]: 


Tf+.x) = - V 2 y,WxVy)S[9iA], (3.1) 

where T is the gravitational potential and C\ parameterizes the strength of the alignment, 
with sign convention such that positive (negative) C\ corresponds to preferential galaxy 
alignment along the stretching (compressing) axis of the tidal held. Positive C\ corresponds to 
the expected alignment behavior and yields an anti-correlation between the intrinsic ellipticity 
of a foreground galaxy and the gravitational shear of a background galaxy. S is a filter that 
smooths huctuations on halo or galactic scales - see section 3.3 for further discussion on the 
effect of smoothing. The gravitational potential is related to the density held through the 
Poisson equation in Fourier space (valid on sub-horizon scales, k cHq): 

Ti A (k, zia) = —4:irGp rn fi(l + zi A )AT 2 5(k, zi A ), (3.2) 

where p rn p is the matter density today. As discussed above, this type of alignment is most 
likely to arise in elliptical galaxies for which angular momentum does not play a significant 
role in determining shape and orientation. 

We can write the tidal alignment model in Fourier space: 

7(+,x)( k > z ) = —Ci(z)p m fi(l + zi A )/(+, X )(k)<5(k, zia), (3.3) 

where we dehne the tidal operator to include the relevant smoothing filter: 

/(+,x)( k ) = sin 2 (9 k ) f { E , B )( k ) S ( k ). (3.4) 

The amplitude C\ will generally be a function of galaxy properties, such as color and lumi¬ 
nosity, as well as observed redshift, rehecting potential IA evolution. The smoothing function 
S(k) corresponds to convolution by a smoothing kernel in configuration space. 
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3.1 Nonlinear contributions 

The tidal alignment ansatz assumes that the intrinsic ellipticity is a purely linear function 
of the tidal field, leaving three sources of nonlinearity which must be simultaneously consid¬ 
ered: (1) nonlinear evolution of the density field; (2) weighting by shape tracer density; (3) 
nonlinear galaxy biasing. These effects are described below. 

3.1.1 Nonlinear density evolution 

We use the methods of standard perturbation theory (SPT) to expand the nonlinear density 
field in terms of the linear density held (see, e.g., [26]). We introduce the notation to 
denote a convolution of n linear density fields with the appropriate gravity kernel capturing 
structure growth (see appendix A.) The nonlinear density held can then be written: 

5 = <5 (1) + <5 (2) + d (3) + • • • (3.5) 

The linear density held is Gaussian, and by Wick’s theorem, correlations of an odd number 
of Gaussian helds are zero. Thus, all correlations are expressed in even powers of the linear 
density held, or equivalently, in an expansion of the linear power spectrum Pi; n . 


3.1.2 Shape tracer density weighting 

Although the tidal alignment model makes a prediction for the (correlated) intrinsic shape 
held as a function of the tidal held (eq. 3.1), the observable quantity that contaminates the 
weak lensing signal is the shape held weighted by the density of the galaxies used as tracers': 
y(x) = (1 + <5 9 (x))7(x). Density weighting arises because we cannot sample the intrinsic 
ellipticity held in an unbiased way, but rather only where luminous tracers form. In Fourier 
space, density weighting results in a convolution (denoted a o b): 7 (k) = 7 (k) + d g (k) o 
7 (k). This weighting has important observational consequences. For instance, although the 
tidal alignment model produces a pure C-mode intrinsic shape held (even when including 
nonlinear evolution), density weighting produces a £>-mode contribution. Density weighting 
also increases measured intrinsic shape correlations on small scales as the contribution of 
galaxy clustering becomes signihcant. 

Finally, density weighting contributes a bias-dependence to the IA amplitude on large 
scales. As seen in appendix A, density weighting produces a term of the form (58/105) 6f (TgPn n , 
where <r| = J P\ in (q)S(q). Due to mode coupling, regions with a high density of shape 

tracers also have large tidal held fluctuations, with an increasing effect for larger tracer bias. 
Density weighting thus yields a bias-dependent enhancement in the observed shape correla¬ 
tions, even if the underlying response to the tidal held is comparable for different tracers. 
This enhancement is present even on large scales (k —> k\[ n ), where linear theory should be 
valid, leading to a bias-dependent effective value of C\\ 


( 58 

PgE{hm) OC Cl f 1 + 

/ CO 

PEE(hin) OC Cf ( 1 T 2 — 


6lCr|J Pun, 

b\(T%\ P] m . 


(3.6) 


It is this effective value of Ci that is actually measured with large-scale intrinsic shape 
correlations. This effect could be largely responsible for the luminosity dependence observed 

7 This situation is completely analogous to observations of the cosmological velocity field. 


-7 - 




in shape correlations (e.g. [5, 6 , 8 , 27]), and the trend with b\ it predicts is consistent with the 
measurements of [ 8 ]. However, this effect may be partly degenerate with other dependence 
of C\ on galaxy properties, including galaxy-halo misalignment [10, 27, 28]. The value of 
(j| depends explicitly on the small-scale cut-off in tidal fluctuations. Moreover, since the 
correction is generally large, higher-loop contributions may be important. Thus, although it 
is possibly to calculate a 2 s under a particular smoothing procedure, it may be preferable to 
instead treat cr| as an additional free parameter of the model. 

3.1.3 Nonlinear galaxy bias 

We include nonlinear bias contributions to the galaxy density field. In Fourier space: 

S g = b\5 + 62(5 o (5) — b2<y 2 + bs(5 o 6 o 5) + ■ ■ ■ . (3-7) 

where <5 is the nonlinear density field, and a 2 = (6 2 ). A more general bias expansion at this 
order would include the nonlocal tidal terms b s and 63^1 (e.g. [29, 30]). We leave consideration 
of these tidal bias contributions for future work. 


3.2 Alignment epoch 

The redshift at which the alignment is set, z\a, and the subsequent IA evolution, are de¬ 
termined by the astrophysical processes involved in galaxy formation and evolution. It is 
sometimes assumed that zia is during matter domination when a halo first forms (“primor¬ 
dial alignment”). However, late-time accretion and mergers could have a significant impact 
on IA, in which case the relevant z\a could be closer to the observed redshift, £ 0 b s - The 
limiting case of z\a = z 0 bs is the “instantaneous alignment” scenario. For further discussion 
on this topic, see [10, 22, 31]. 

At leading order, the impact of alignment epoch can be absorbed into a redshift- 
dependent amplitude, A(z), and eq. 3.3 can be written: 

7 (+, x )( k > z ) = A (z)f(+,x)( k)<S(k,z). (3.8) 

In the case of instantaneous alignment (21 a = z 0 bs): 

A(z) = -Ci(z)pm,oO- + z). (3.9) 


For earlier epochs of alignment, the linear order evolution of the density field evolution can 
be absorbed into the growth factor D(z ): 

A(z) = -Ci{z)p m , 0 (1 + ,ia)^. (3.10) 

In the case of primordial alignment during matter domination, the IA held is insensitive (at 
linear order) to zja- 


A(z) 


-Ci(z)p mfi (l + z) 
D(z) 


(3.11) 


where D = (1 + z)D(z) is the rescaled growth factor, normalized such that D(z) = 1 during 
matter domination. In these expressions, we have separated redshift dependence due to the 
growth of structure from possibly evolution in the response to the tidal held, captured in 
Ci(z). 
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Once nonlinear evolution of the density field is included, eq. 3.8 is only exact in the case 
of instantaneous alignment. Equations 3.10 and 3.11 provide only approximate amplitude 
scalings when z\a > ^obs, since the choice of z\a will affect both the overall amplitude of 
the IA field as well as its scale dependence (via nonlinear evolution). In this case, nonlinear 
predictions for the observed IA correlations must be treated carefully. For instance, if the 
intrinsic shapes and alignments of galaxies are set at the redshift of initial formation, zp. 
the tidal field should be evaluated at that redshift. A model is then required to describe 
the subsequent evolution of halo and galaxy alignments. However, the relevant shape and 
density tracers fields should still be evaluated at £ 0 b s - 

Moreover, even if the alignment field is set at zp and simply “passively” evolves to 
z 0 bs, it will be effectively transformed by the motion of the halos/galaxies on which it was 
initially imprinted. In our Eulerian approach, the advection of modes due to fluctuations 
on larger scales will suppress correlations of density (or tidal) fields at different redshifts, 
which is non-physical in the case where primordial alignments are “carried” with the moving 
galaxies. At one loop order, this effect appears as a relative suppression of the mode-mixing 
term compared to the propagator contribution. 

One simple approach to modeling the evolution of IA set at high z\a has been to take 
the geometric mean between the power spectra evaluated at the two relevant redshifts, zja 
and £ 0 bs (e.g. [22]). While such an approach avoids the non-physical suppression that arises 
in the Eulerian treatment, it also fails to fully treat advection effects and should thus be 
treated with caution. In this work, we adopt the instantaneous alignment ansatz, for which 
Eulerian SPT is well-suited. We leave a more thorough treatment of nonlinear advection of 
the IA field for future work. 

3.3 Smoothing the tidal field 

Tidal field fluctuations on sufficiently small scales should not have a significant impact on 
the observed IA of halos/galaxies, and it is thus common to smooth the tidal field (e.g. 
[9, 15]). The correct prescription for smoothing is not settled - from a physical standpoint, 
smoothing reflects the finite size of the halos/galaxies and the regions from which they formed. 
Alternatively, the smoothing can be viewed as an element in any effective theory, in which 
small-scale physics is absorbed into the measured parameter values. 

We treat the smoothing filter as a multiplicative function in Fourier space. Smoothing is 
applied to the tidal field, via a redefinition of the IA operator (eq. 3.4), before the convolution 
integrals are performed. We choose the form: 

S(k) = e (^m) ; 

, 2vr 

k S m = a—, (3.12) 

Kl 

where Rp is the Lagrangian radius of the halo, and a is an order-unity factor to capture the 
uncertainty in the IA response to smoothing. The Lagrangian radius of a halo with mass 
Afhaio is defined as the radius of a sphere in a homogenous universe that contains the halo’s 
mass: 

47T Q 

ALhalo = PmR 3 L■ (3.13) 

In the tidal alignment picture, it is often assumed that halo shapes are determined by the 
external tidal field during collapse, making the Lagrangian radius a natural scale for smooth¬ 
ing. Fluctuations on smaller scales may significantly influence the actual shapes of the inner 
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regions of halos (and the galaxies residing there), but should not contribute to correlations 
on large scales [10]. For the typical luminous red galaxies (LRGs) used for comparison in 
this work, the halo mass is Mh a i 0 ~ 3 x 10 13 h^ 1 M & . We can write Rl, scaled by typical 
parameter values: 

R L - (4.41 ft-Mpc) x (_J^L _) 1/3 (^)- 1/3 . ,3.44) 

Note that in comoving coordinates, Rl is constant with redshift, since the product of volume 
and density remains fixed. In this work, we choose k sm = 1.0 h Mpc -1 . The mass dependence 
of the Lagrangian radius, and this smoothing scale, is fairly weak. 

The tidal field on smaller scales may have an important role in determining actual galaxy 
shapes and orientations, particularly when considering the central baryonic component, but 
these effects simply modulate the large-scale correlations, analogous to how nonlinear galaxy 
formation physics modulates the large-scale galaxy bias. The physics of IA may actually 
yield a different smoothing scale than the one assumed here. Indeed, comparing predictions 
with different smoothing scales to measured IA provides one way to probe the underlying 
effects of tidal field fluctuations on different scales. 

3.4 Redshift-space distortions 

The peculiar velocities of galaxies lead to shifts in their apparent distance along the line-of- 
sight. On linear scales, the effect of these “redshift-space distortions” (RSDs) on the galaxy 
density field are well-described by the Kaiser formula [32]: 

S s g (k) = (b l + f f , 2 )S r (k), (3.15) 

where / = din D/ din a is the logarithmic growth rate, n = k ■ z, and “r” and “s” denote 
quantities in real and redshift space, respectively. Although the apparent line-of-sight posi¬ 
tions of galaxies are affected by RSD, their observed shapes are not. Thus, at leading order: 
7 s = 7 r - The density-weighted intrinsic shape field in redshift space ( 7 s ) is affected by the 
shift in observed positions of shape tracers as well as higher-order effects on r y s . For further 
discussion, see the appendix of [ 8 ]. 

For the projected statistics considered here, only RSDs on linear scales will contribute 
significantly, suggesting that we can consider only the linear contribution to 7 s , for which 
the effect of RSDs can be ignored. In this case: 

PgE = PgEi 

Pee,bb = Pee,bb■ (3.16) 

We use this correction in the results presented below. Note that the treatment of RSDs in 
[10, 15] assigns a linear-order effect to both 5 g and 7 , yielding a small misestimate of the IA 
correlation function. For instance, for a projection length of II max = 80/i _1 Mpc, the RSD 
treatment of eq. 3.16 yields to a ~ 2% correction to the real-space w g + at r p = 50/i^ 1 Mpc, 
compared to a ~ 4% correction in [10, 15]. At smaller projected separation, RSDs rapidly 
become negligible. 

The linear theory result of eq, 3.16 neglects an 0(Pj? n ) contribution from density weight¬ 
ing that yields a term proportional to P] ni on large scales (see section 3.1.2). This term is 
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present in both P^ E and P EE and will contain an additional Kaiser-like enhancement factor. 
However, given the small impact of RSDs on the scales of interest, we do not consider this 
contribution here. 


3.5 Combining all contributions 

In this section, we present the key results from combining all relevant contributions to tidal 
alignment. Details can be found in appendix A. Using SPT at one loop, we calculate the 
relevant power spectra, defined in eq. 2.5. We then transform to the correlation functions 
(eq. 2.8) to compare with measurements. 

Combining eqs. 3.7-3.8 and keeping terms up to 0(6 4 ), the density weighted intrinsic 
shape field is: 

7 i(k, z) = A(z) ((1 - b- 2 cr 2 )fi6 + bid o (fid) + b 2 (S o 5) o (fiS) + b 3 (d o 6 o 6) o (fid)) (3.17) 

Collecting all contributions to the cross-correlation between galaxy density and shape that 
are not zero by symmetry, we find: 


(«5 9 |7+) = A(z) x [hf+m + bl(d\d o (f + 5)) + b 2 f + (d o 5\S) 
+ b 1 b 2 (dod\do(f + d)) c + 0(6 6 L )], 


(3.18) 


where the c subscript indicates connected correlators, with two or more branches between 
the two points. Similarly, the intrinsic shape auto-correlations can be calculated, up to one 
loop: 


<7i|7i> = A 2 (z) x f 2 (d\d) +2^(616 o(f l 5)) + b 2 (do(fid)\do(f i d)) c + 0(dl) 


(3.19) 


We define power spectra corresponding to the correlators in equations 3.18-3.19, recalling 
that we can evaluate these correlations at k y = 0: 


P gE (k) — A(z) sin 2 (0k) [&i-fo|S + b 2 Po\oE + b 2 P 00 \s + bib 2 P 00 \ 0E + 0(P E )\ , 
Pee( k) = A 2 (z) x [sin 4 ((9 fc )Ps|s + 2h sm\6 k )P S \ 0E + b 2 P OE]OE (0 k ) + 0(Pl)\, 
Pbb( k) = Al 2 (z) x \b 2 P 0Bl0B (6 k ) + 0(Pl) 


(3.20) 

(3.21) 

(3.22) 


The subscripts in these expressions indicate the fields being correlated: density field (0), 
smoothed density field (S), or components of the smoothed tidal field (E or B). For instance, 
Pq|5 is the nonlinear matter power spectrum in which one power of the density field is 
smoothed at the relevant IA scale. Dependence on 0 k is explicitly indicated. For the exact 
form of these terms in SPT, see appendix A. When the projection length is larger than 
the scale of interest, the Limber approximation can be used, and only k z = 0 modes will 
contribute, corresponding to sin 2 ( 14 ) = 1. 

Although we arrived at equations 3.20-3.22 through an SPT expansion, each term can 
be evaluated to arbitrary precision. For instance, SPT is known to provide a comparatively 
poor fit to the nonlinear matter power spectrum, (<j|<5) (e.g. [33]). More accurate techniques 
are available to calculate this term (e.g. a fitting formula such as Halofit [34] or a simulation 
emulator [35]). In the following results, we use a Halofit nonlinear power spectrum for P 0 |0i 
P 0 \ S , and P s{s . 
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Figure 1 shows the linear and nonlinear contributions to w g + and w++ in the tidal 
alignment model for an LRG-like shape and density tracer sample (e.g. [28, 36]), for which 
we assume bias values b\ = 2.1 and 62 = 0.5. The Ogi^m contribution from density weighting 
rescales the linear theory prediction (i.e. the effective value of C\ at large scales) - it is shown 
with an amplitude set by the one-loop prediction (see section 3.1.2). While the contributions 
from nonlinear galaxy bias are subdominant, the effect of density weighting is significant and 
is larger than the correction from nonlinear dark matter clustering for these highly biased 
tracers. We thus confirm that the NLA model omits important nonlinear effects. 

3.6 Comparison to observed IA 

We compare predictions of this nonlinear tidal alignment model with two different IA mea¬ 
surements using highly biased tracers with spectroscopic redshifts from the Sloan Digital 
Sky Survey (SDSS, [37]). The measurements of [28, 36] use ~ 84,000 luminous red galaxies 
(LRGs) from SDSS DR 6 , with spectroscopic redshifts in the range 0.16 < z < 0.47 and mean 
redshift z = 0.32. We also compare with the recent measurements of [ 8 ] from the BOSS 
“LOWZ” sample from SDSS DR11. Within the redshift range 0.16 < z < 0.36 (for which the 
sample is roughly volume limited), [ 8 ] use ~ 160, 000 galaxies with good shape measurements, 
yielding a mean redshift z = 0.28. 

While the galaxy properties in these two samples are similar, the LRG sample used by 
[28, 36] includes objects at higher redshift and has a higher average bias than the LOWZ 
sample of [ 8 ]. In addition, the LRG measurements of [28, 36] use isophotal galaxy shapes 
which have not been corrected for the point-spread function (PSF). Isophotal shapes apply 
more weight to the outer regions of galaxies, leading to larger measured values of ellipticity 
and all related statistics . 8 For the IA autocorrelation, the LRG measurements of [36] 

were reprocessed by [ 10 ] to account for projection effects and the treatment of randoms in 
the original estimator. The LOWZ w++ measurement of [ 8 ] has significantly lower signal-to- 
noise. This discrepancy may be due to the use of different shape measurement techniques or 
the projection procedure used in [10]. See [ 8 ] for further discussion. 

Figures 2-3 compare the tidal alignment model with all 0(Py m ) corrections with the 
shape correlation measurements of these two different samples. The alignment amplitude is 
fit to measurements of w g + above 2 h _1 Mpc are shown for both fixed erf, (calculated using 
by smoothing the density held as discussed in section 3.3) as well as allowing cr| to be a free 
parameter (constrained to be positive). Also shown are the linear and NLA model predictions, 
with and without smoothing of the tidal held, ht to measurements above 10 /? _ 1 Mpc. The 
model is separately ht to the LRG w ++ . As discussed in [10], the amplitude on large scales 
is consistent with that of w g ±. The low signal-to-noise in the LOWZ w++ measurement 
prevents a meaningful model ht. For reference, the linear and NLA models, normalized using 
the ht to the LOWZ w g +, are shown in hg. 3. In summary, including all nonlinear corrections 
improves the ht when compared with the linear or NLA models. Allowing cx| to be a free 
parameter further improves the ht, particularly for the LRG w++ measurements. 

4 Tidal alignment on small scales 

The perturbative treatment described in the previous section provides a consistent approach 
to understanding the contributions to shape correlations in the tidal alignment picture. How- 

S R. Mandelbaum, private communication. It remains an open question whether the PSF introduces spu¬ 
rious correlations that contaminate w g + and w++. 
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Figure 1. The linear theory result, NLA correction, and all additional 0{P'y m ) corrections are shown 
for w g+ (top panel ) and ui++ (bottom panel), with by = 2.1 and &2 = 0.5. Overall normalization 
for both statistics, determined by Cy, is chosen to match the LOWZ w g+ measurement of [8], while 
the value of er| is calculated using the smoothing discussed in section 3.3, with fc sm = 1.0 /iMpc’ 1 . 
Contributions are labeled using the associated pre-factors in eqs. 3.20-3.22, and the <JgPn n contribution 
(black) is separated from the other tracer density weighting terms (green). Note that the “NLA corr.” 
term shows the difference between the NLA (Halofit) and linear theory predictions. Dashed lines 
indicate a negative value. 


ever, it will necessarily break down on sufficiently small scales. In particular, the tracer 
clustering (density auto-correlation) on small scales is not well-described in a perturbative 
framework. In this section, we employ an alternative approach, explicitly separating the 
effects of galaxy clustering and galaxy alignment and modeling each separately. 

Without loss of generality (assuming approximately constant source number density 
across the projection length, 2II max ), we can divide the shape correlation statistics into 
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Figure 2. Data points show the LRG measurements of w g + (top panel) and w++ (bottom panel) from 
[28, 36]. The solid lines show the best fit linear, NLA, and total model predictions, with bi = 2.1 and 
b 2 = 0.5. For comparison, dashed lines indicate the linear and NLA models with no smoothing of the 
tidal field. The alignment amplitude and cr| is fit separately to w g + and w++. 


contributions from the intrinsic shape held and galaxy clustering: 

W a+( r p) = ( r r+) t r p aC ( w 99( r p) + 2II max), 

w (++,xx)( r p) = (7(+,x)7(+,x))r p ( w gg( r p) + 2II max ). (4.1) 

Equation 4.1 defines the quantities ( 7 + )^ ac and ( 7 ( +iX ) 7 ( +iX ))^ ac , which are interpreted as 
the average alignment between shape and density tracers and between two shape tracers, 
respectively, at a given projected separation. As this division makes clear, on large scales, 
the clustering of shape tracers can be ignored ( w gg <C 2II max ), and shape correlations are 
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Figure 3. Data points show the LOWZ measurements of w g + (top panel) and w++ (bottom panel) 
from [8]. For w g +, the solid lines show the best fit linear, NLA, and total model predictions, with 
&i = 1.8 and & 2 = 0.5. For comparison, dashed lines indicate the linear and NLA models with no 
smoothing of the tidal field. The linear and NLA models for w++ are shown, normalized using the fit 
to w g + due to low signal-to-noise in w ++. 


determined entirely by these average quantities. On smaller scales, the contribution from 
tracer clustering can be significant. 

4.1 Modeling ( 7 + )^ ac around halos using tidal alignment 

In this section, we construct a model for ( 7 +}j? ac around biased tracers of dark matter (e.g. the 
correlations relevant to the w g+ statistic). Similar arguments apply to modeling rc( ++iXX ). 
The average of the intrinsic shape field at projected separation r p from a density tracer can 
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be written: 


( 7+ ) , p = ^— [ dU(6 g \ 1+ ) r (4.2) 

n ma x 

Note that this expression is very similar to the definition of w' 5 +, except the effect of weighting 
by the shape tracer density is not present. 

In the tidal alignment framework, the intrinsic shape field in a given location is de¬ 
termined by the local tidal field. The quantity ( 7 +) is thus determined by the correlation 
between density tracer positions and the tidal field. Following the results of section 2, it 
is clear that ( 7 _|_) will essentially be the projected cross-correlation between galaxies (the 
density tracer) and dark matter, except that the derivatives relating the density to the tidal 
field will result in a J 2 rather than a Jo Hankel transform: 

(7+) = — A(z)w gS , 2 , 

00 

1 A(z) f 

~ -777TV / dKKJ 2 {nr p )P g s(K), 

Climax (^TT) J 

0 

where the second line assumes the Limber approximation (e.g. eq. 2.9), and we have intro¬ 
duced the correlation function notation: 


(4.3) 


Cab,n{ r p ) II) 


1 


00 

J dk z d,KK cos(k z U)J n (Kr p )P ab (k, 9 k ) 
0 


(4.4) 


for the correlation between quantities a and b, weighted by the relevant Bessel function J n 
to account for the azimuthal angular dependence (see eq. 2.7). 

To relate ( 7 +) to ( 7 + ))l ac we must include the effect of sampling the tidal field at the 
locations where shape tracers form. As discussed in section 3.1.2, in the tidal alignment 
framework, this effect produces an additional contribution to alignment ~ b\o 2 s P\\ n , which 
changes the effective large-scale alignment amplitude: = (1 + 58/1056ifj|)C'i. We will 

assume that ( 7 + ))l ac includes this additional term, such that on large scales ( 7 _|_}* 4ac = 
(Cj ;ff /C'i)( 7 +} ? , p . The contribution from shape tracer clustering has been explicitly separated 
into the w gg term in eq. 4.1. 

<7 + >r=<7 (4 - 5) 

We are free to choose a model for P g $, the cross-power spectrum of galaxies and dark 
matter. This correlation could be treated perturbatively, as done in the previous section, 
by including nonlinear galaxy bias. Instead, to describe tidal alignment down to highly 
non-perturbative scales, we use a Navarro-Frenk-White (NFW) halo density profile [38]: 

P(r) = (r/r.)(l + r/r.r ^ 

defined by its scale density ( p s ) and radius (r s ). Equivalently, the NFW profile can be 
described with the halo mass (typically defined at some density threshold) and concentration 
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(the ratio of the profile scale radius to halo radius, defined by the same density threshold). 
Using galaxy-galaxy lensing, and assuming the mass-concentration relationship of [39], [8] 
measure the average LOWZ halo mass to be Migoft = 1.51 x 10 13 /j - 1 M q , where Mi§o 6 denotes 
the mass enclosed within a region of 180 times the average matter density. To model P g s, we 
truncate the NFW profile at 2Ri80b, take its Fourier transform, and combine with the linear 
P g s(k) = biP\- m (k) to include correlations on larger scales. We note that more sophisticated 
models for P g s could be used within this framework, particularly to describe the transition 
region between the one- and two-halo regimes. Nevertheless, our approach provides relatively 
accurate predictions. Rather than apply a smoothing filter to the tidal field, we explicitly 
separate the alignment behavior within the one-halo regime (discussed below in section 4.2). 

Modeling w g+ requires a description of both (7 + )]4 ac and w gg . An analytic description 
of Wgg on small scales, particularly the one-halo regime, requires halo modeling specific to 
the relevant shape and density tracer samples. However, we are also free to use an empirical 
model (from simulations or observations). In general, the projected clustering signal will be 
measured to significantly higher precision than the associated shape correlations, and thus 
in some contexts, using the observed w gg as a component of the intrinsic alignment model is 
feasible. On small scales, w gg 3> 2n max , and thus w g + ~ ('y + ) t ™ c w gg . Within the one-halo 
regime, as we argue below, (7+)[I ac may be approximately constant, leading to w g+ oc w gg . 

Figure 4 compares the result of our model for ( 7 + ))f ac with the LOWZ observational 
results from [8]. With the treatment of the one-halo regime discussed in section 4.2, our 
model for ( 7 + )[f ac agrees well with the observational results. For comparison, we also show 
the SPT results of section 3. In the case where density weighting is included, we calculate 
(7 + )]4 ac by dividing the model prediction for w g+ by (w gg + 2n max ), where we have used 
Halofit nonlinear matter clustering and linear bias to model w gg . For the NLA model, which 
does not account for source density weighting, we show w g+ divided by 2n max only. 

4.2 The one-halo regime 

On large scales, in the two-halo regime, alignments largely reflect initial correlations in the 
tidal field. On these scales, tidal alignment (and perhaps tidal torquing) mechanisms should 
accurately describe intrinsic shape alignments. In the one-halo regime, however, information 
of these initial conditions is largely scrambled by nonlinear gravitational processes. We thus 
expect that these theories will no longer provide a qualitatively correct description of shape 
alignments, even if a fully non-linear treatment of the tidal field (e.g. an NFW halo profile) 
is used. On these scales, ( 7 +))f ac should differ from the tidal alignment prediction oc w g s t 2 
made above. 

Understanding the alignment of central and satellite objects within the same halo re¬ 
mains a crucial question in both observational and theoretical studies. Within a halo, there 
are three potential types of shape-position alignment: (1) alignment of central shapes with 
satellite positions; (2) alignment of satellite shapes with the central position; (3) alignment of 
satellite shapes with other satellite positions. The alignment of central shapes with satellite 
positions has been measured (e.g. [8, 40]). Since satellites should essentially trace the halo 
density, this effect is expected in the tidal alignment picture and may be related to recent 
observations showing that satellite galaxies are preferentially located along the axis of sur¬ 
rounding filaments [41]. Conversely, satellite alignment appears to be significantly weaker 
[42], although the most luminous satellites may be somewhat aligned with central positions 
and the surrounding large-scale structure [8, 43]. 


17 - 



Figure 4. Model results for ( 7 + )j. rac , shown for an NFW profile for M 180 {, = 1.51 x 10 13 h 1 M 0 , is 
compared with the observational results from the LOWZ sample [8]. The dashed blue line shows the 
continuation of the tidal alignment prediction below R^ooc (indicated by the vertical line). Predictions 
from the full SPT model of section 3, as well as the NLA model (dashed line indicates no smoothing), 
are shown for comparison. As discussed in the text, the NLA model does not account for density 
weighting, and the associated prediction for ( 7 + )jf ac does not include w gg . 


A detailed description of these intra-halo alignments is beyond the scope of this paper. 
Instead, motivated by the observational results of [8] (see figure 4), we will make the following 
simple ansatz: within the halo ( 7 +)^ ac saturates at a constant value. The radius at which this 
saturation occurs will in general depend on the properties of the shape and density tracers 
under consideration. It can be related to the halo virial radius or be left as a free parameter. 
In this work, we set ( 7 _|_)*^ ac to a constant value at i? 200 ei defined by mean halo density 
p(r < i? 200 c) = 200/? cr it - This value is somewhat inside the halo virial radius and corresponds 
to R 200 C = 0.43/i _1 Mpc for the average LOWZ galaxy halo. Note that the effect of projected 
correlations from pairs of galaxies not in the same halo rapidly drops for r p < -Rhaioj and 
thus the pairs of objects that dominate the measured correlations on these scales are indeed 
in the same halo. 

If ( 7 + )^ ac is indeed dominated by central shape alignment with satellite positions, sat¬ 
urating at a constant value implies that the average satellite distribution maintains a rea¬ 
sonably constant shape and orientation as a function of radius. The apparent lack of scale 
dependence within the halo may be in tension with the behavior of dark matter in simu¬ 
lations, which exhibits significant internal “twisting” (e.g. [44]). One possibility is that the 
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satellite distribution may maintain a more regular shape than the dark matter. Alternatively, 
the behavior of ( 7 _|_))f ac seen in the LOWZ galaxy sample may not be typical of less massive 
satellites. Potential systematic uncertainties due to fiber collisions and shape measurement 
could also affect (7+))4 ac measured on these small scales (see [8]), in which case the observed 
saturation may not be physical. In any case, some degree of saturation or other departure 
from tidal alignment is necessary to avoid divergent alignment on small scales. 

The small-scale alignment behavior of LOWZ galaxies may not be universally valid when 
considering different halo masses and satellite properties. Further understanding of shape 
alignments within halos remains an important goal. We suggest that the ( 7 _|_))f ac statistic, 
as well as similar statistics such as those used in [8, 45], be used to probe alignment around 
centrals and satellites with different masses and other properties. For instance, the results 
presented in figures 11-12 of [45] show that average alignments of galaxy shapes resolved in 
hydrodynamical simulations reach a similar saturation within the one halo regime. Further 
observations ands simulations, including centrals and satellites spanning a wide range in color 
and luminosity, will be helpful in probing these questions. 

4.3 Contribution to cosmic shear 

IA contributes a potentially important systematic signal to the measured cosmic shear. In 
previous sections, we have considered IA correlations between biased tracers. However, the 
relevant intrinsic alignment contamination to cosmic shear measurements is sourced by corre¬ 
lations between the shapes of tracers and the dark matter density field itself. In this section, 
we will briefly describe tidal alignment modeling of these terms. For a more detailed study 
of the impact of IA contamination for various lensing surveys, see [46]. 

The two-point shape correlations can be expressed in terms of the E and B modes of 
the angular power spectrum (or angular correlation function): 

C EE = c EE,GG + C EE,GI + C EE,n 

C BB = C BB ’ n , (4.7) 

where the GG term is the standard cosmic shear signal, GI is the cross-correlation between 
lensing shear and IA (which remains even when cross-correlating shapes from different redshift 
bins), and II is the IA auto-correlation from pairs of galaxies that are affected by the same 
fluctuations in the gravitational potential. The GI term corresponds to (<5|7_|_), while the II 
term comes from (7(+ )X )|7(+,x))- See [9] for further discussion. 

In our perturbative (SPT) treatment of tidal alignment, these correlations are essentially 
the same as those presented above (section 3.5), with the galaxy density field 5 g on the left 
side of the GI correlations replaced with the dark matter density <5 (the 5 g appearing in the 
tracer density weighting remains). Comparing to eq. 3.20, we find at one loop: 

Pse( k) = A(z) sm 2 (8 k ) [P 0|S + hP 0l0E + 0(Pl)\ . (4.8) 

The II correlations of eqs. 3.21-3.22 are unchanged. See appendix A for the derivation of 
these terms. 

For our non-perturbative treatment, explicitly separating the effects of alignment and 
clustering, we must now consider the average alignment of shape tracers sourced by the dark 
matter density field itself. While the relevant IA auto-correlation (II term) is unchanged 
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from equation 4.1, the IA-density cross-correlation (GI term) in configuration space is: 

^<5+ = <7+)r p M ( 2n max + W gS ), 

= ^~ A ( Z ) + Y^ b ' a s wl sl?j (2n max + w gS ). (4.9) 

The quantity ( 7 + )^ M is the average IA correlation sourced by dark matter density, measured 
at the positions of shape tracers, which leads to the same large-scale enhancement found in 
eq. 4.5. Although not directly observable in galaxy shape measurements, this ws+ correlation 
is closely related to the GI contribution to cosmic shear and can be directly measured in 
simulations (e.g. [45]). 

Figure 5 shows the tidal alignment predictions for w$+ (the GI contribution) from 
shape tracers similar to the LOWZ sample for both the SPT model and the non-perturbative 
(NFW halo model) approach. For the SPT modeling, we continue to smooth the tidal field at 
k sm = 1 /iMpc -1 . In the non-perturbative model, the tidal field is not explicitly smoothed, 
although following the treatment in section 4.2, the value of ( 7 + ))f ac is set to a constant 
within .Rconst = -R 200 c- The prediction for w g s is made using both Halofit (with linear bias) 
and an NFW profile. For illustrative purposes, we consider only dark matter halos at the 
LOWZ mass. This approach is overly simplified - further investigation on alignments within 
the one-halo regime as well as including contributions from the entire halo mass range are 
required. 

5 Discussion 

In this work, we have developed the tidal alignment framework for galaxy intrinsic align¬ 
ments. We included nonlinear contributions from source density weighting, dark matter 
clustering, and biasing. We have also explored the effect of smoothing the tidal field on the 
IA correlations. Accounting for these effects improves the consistency of these models beyond 
the linear and nonlinear alignment (NLA) models and provides a more accurate description 
of the observed alignments of luminous red galaxies. In their fiducial forms, these models are 
constructed without introducing additional free parameters (beyond the IA response value 
C\ and the galaxy bias parameters or relevant halo mass). The other required elements - the 
smoothing scale and resulting value of cr 2 s and, in the case of the one-halo term, the radius 
where ( 7 +) becomes constant - can be predicted from physical arguments or left as free pa¬ 
rameters within a narrowly allowed range. Since marginalizing over an increased number of 
free parameters in the IA model will reduce the constraining power of lensing measurements, 
minimizing additional free parameters is an important advantage to our models. 

Of the three sources of nonlinearity that contribute in the tidal alignment picture, source 
density weighting is the most significant, particularly for highly biased tracers. Density 
weighting effects arise because we observe intrinsic alignments only at the locations where 
galaxies form. In the case of tidal alignment, we are thus sampling the tidal field in these 
special locations. The clustering of shape tracers increases IA correlations on small scales and 
generates a .B-mode autocorrelation signal. Somewhat less intuitively, the impact of density 
weighting extends to large scales, where it effectively enhances the amplitude of the linear 
theory prediction. This term contributes to an effective value of C\. which depends on both 
the bias of the shape tracer and the relevant smoothing scale of the tidal field. Physically, 
this effect is a direct result of sampling the tidal field in the special locations where galaxies 
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Figure 5. Predictions for w$+ are shown for both the SPT and non-perturbative models. Shape 
tracers are assumed to have the same bias and alignment properties as the LOWZ galaxies [8]. As 
discussed in the text, (7+) saturates below r p = -R2000 and results are shown for w g $ predicted using 
both Halofit (with linear bias) and an NFW profile. For reference, the (7+) line (blue) shows the 
alignment effect without shape tracer clustering. 


form. More highly biased objects tend to form in regions with larger tidal field fluctuations, 
leading to larger intrinsic ellipticity correlations, even if the underlying response to a given 
tidal field is the same for all halos/galaxies. 

We have also shown that observed shape correlations can be modeled non-perturbatively 
by separately considering lensing and alignment components. In particular, we developed a 
straightforward description of alignment on small scales using the tidal field from an NFW 
halo profile and saturating the value of ( 7 +)^ ac inside the halo. This model is likely overly 
simplistic - further work in modeling, observations, and simulations are required to better 
understand these effects. Nevertheless, this model is simple to implement and likely to be 
more accurate than existing models on smaller scales where a significant amount of lensing 
information is contained. 

The two approaches described in this work are complementary and may be useful in dif¬ 
ferent regimes. The SPT model is straightforward to implement and provides a consistent and 
well-understood description on intermediate and large scales. Splitting the correlations into 
alignment and clustering components provides more flexibility in including non-perturbative 
physics, which allows a description on small scales, including within the one-halo regime. 
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The model’s accuracy is then determined by the particular choices for the nonlinear tidal 
field and galaxy/matter clustering. Both of these approaches can be generalized to include 
other alignment mechanisms (such as tidal torquing). 

As plans are made for analyzing current and future weak lensing experiments, it is im¬ 
portant to emphasize that the simplest IA models may not correctly describe all galaxies. 
Assuming tidal alignment when the true astrophysics of galaxy formation produces more 
complicated behavior will yield systematic bias in the final parameter constraints. Moreover, 
even forecasting the impact of IA requires a model both sufficiently accurate and complex 
to avoid overestimating the statistical power of a given survey and mitigation scheme [46]. 
In this paper, we have neglected contributions from higher order dependence on the tidal 
field (i.e. tidal torquing models). Although symmetry allows such terms, and they are likely 
present at some level for galaxies kinematically dominated by angular momentum, current 
observations provide no strong evidence that these terms are significant. Nevertheless, until 
higher-order alignment contributions can be quantified or more tightly constrained in obser¬ 
vation or realistic simulation, it will be important to include them. In future work, we plan 
to integrate these terms into a consistent framework. 

The tidal alignment picture has been successful in describing IA in both astrophysical 
observations and simulations. In consistently including nonlinear contributions, this work 
should improve the accuracy and flexibility of tidal alignment as a tool in precision cosmol¬ 
ogy. Furthermore, studying these contributions probes astrophysical effects which must be 
understood as we develop a more complete model of intrinsic galaxy shapes. 
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A Calculating nonlinear contributions in SPT 


Details of the calculations for P g E, Pee, Pbb are presented below. 


A.l Calculating P g E 

We start by determining all possible correlations which can contribute to P g E- 


(<y7+) =A(z) x 


(1 


b 2 a 2 )b l f + {5\5) + b 2 (5\6 o ( f+5 ])) + b 2 f+(5 o 5\5) 


+ hhiSKS o 5) o (/+<*)> + b 3 f+((6 o <5 o <5)|<y) + hb 2 (S o 5\5 o (f+S)) + ... 


(A.l) 


where dots represent the terms that contribute at higher than one-loop order. We remove 
terms that are zero by symmetry, and by decomposing the remaining correlators we find: 


(^317+) = A ( Z ) x 


((1 - b 2 a 2 )b 1 + b±b 2 a 2 + 3b 3 a 2 )f+(5\6) + bj(S\5 o (f+5)) + b 2 f+(5 o <f|5> 


+ hhmS o 6) o (f+S)) c + b 3 f+((5 o <5 o 5)|(5) c + hb 2 (S o 5\5 o (f+5)) c + ... . 

(A. 2 ) 


The “c” subscript denotes the tightly coupled component of a correlator (i.e. components 
not proportional to a 2 ). After the bias renormalization b\ —>• b\ + 3b 3 a 2 (e.g. [47, 48]), and 
keeping all remaining one-loop terms (i.e. 0(6f in )) we find: 


(^ 3 17+) = A(z) x hf+m + b 2 (5\6 o (f+5)) + b 2 f+(d o 5\5) 

+ b 1 b 2 (5o5\5o(f+5)) c + 0(5& l j ■ 


(A.3) 


Although derived through perturbative considerations, this model for P g E can now be consid¬ 
ered more generally, and each term can be evaluated to arbitrary precision. This is possible 
since each term has a unique bias dependence, and the corresponding correlators are IR safe 
[49]. While the particular bias renormalization scheme will change depending on the choice 
of PT technique (or loop order), the overall form will not. 

We now consider the correlators from equation A.3. Note that the resulting convolution 
integrals often include the IA operators /+ and / x . In appendix B, we derive the relevant 
angular integrations over these operators. 


A.1.1 The (<5|<5) term 
Expanding this term at one-loop order: 

(£|5) = (d (1) |5 (1) > + (5 (2) |d (2) ) + 2(<5 (1 V (3) }- (A. 4) 

In terms of the power spectrum with the standard definition, (27r) 3 P(A:)o|o(I i: '(k-|-k / ) = (<5|<5), 
we have: 

Po\o(k) = P L (k) + 2ioo(fc) + 6/c 2 J 00 (k)P L (k), (A.5) 

where Too and Joo are the one-loop SPT corrections to the matter power spectrum (e.g. [50]): 

/»(*) = / ^3 P L(9)-PL(|k-q|)F 2 2 (q,k-q), (A.6) 

1 f d 3 a 

J oo( k ) = -p J ^3^L(9)^3(k,q,-q). 
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We have employed the relevant gravity kernels: 


F 2 (q,k —q) 
F 3 (k,q, -q) 


7// + 3 a — 10a// 2 


14a(l + a 2 


3024 q 2 


12 


a^ 


2 a//) ’ 

158 + 100a 2 



(A.7) 

i)S < 7 “ 2+2 M^)/ 


where a = q/k and // = k ■ q. The smoothing of the density field easily factors out, yielding: 
P 0 \ s (k) = S(k)P Q \ Q (k) and P s \g(k) = S 2 (k)P 0 \ 0 (k). Although this expression for P 0 |o(^) is 
correct at one-loop order, it is known to provide a poor description on comparatively large 
scales. In this work, we instead use the Halofit fitting formula [34] for nonlinear P 0 |o(^)- 


A.1.2 The (<5|d o ()) term 
Expanding the term to one-loop order: 

(S\S o ( f + 5 )) = <<5 (2 V (1) o (f+sW)) + <£ (1 V {2) ° (,M {1) )> + (<5 (1 V (1) o (/+<5 (2) )>. (A.8) 


These correlators correspond to contributions to the power spectrum: 

Po\o E (k) = A 0 \ 0E (k) + B 0 \ 0E (k) + C 0 \ 0E (k). (A.9) 

These terms can be written: 

A )|0 E {k) = 2 J -^^f + (-q)F 2 (k-q,q)P L (q)P L (\k-q\), 

B 0]0E (k) = 2P L (k) J ( 03/ + (-q)P 2 (-k,q)P L (g), 

C 0 \ 0E (k) = 2 P L {k) J ^f + (-k + q)F 2 (-k,q)P L (q). (A.10) 


Using eq. B.2 for the integral over /+ and preforming the coordinate transformation q —> k—q 
in C Q \ 0E {k) yields: 

Ao\oe(^) = 2sin 2 (6 ) fc) J ^-^P 2 (/j, q )S(q)F 2 (k - q, q)P L (q)P L (\k - q|), 

Bo\OE(k) = sm 2 (9 k )^cr 2 s P L (k), 

C 0 \ 0E (k) = 2 S hi 2 (9 k )P L (k) J -^V 2 (fi q )S(q)F 2 (-k,k-q)P L (\k-q\), (A.ll) 


where a 2 = 

P(.k)o\ 0 E 


/ (C)3 PUq) and a 2 = J 


d 3 q 


58 


v 2 S P L(k) + 2 


105 
+ 2 P L (k) 


d 3 q 

(2ir) 3 


(2tt) 

S(q) 


d 3 q 

(27t) 3 


PL(q)S(q). Combining these terms, we have: 


zP 2 {Hq)S{q)F 2 (k - q, q)P L (q)P L (\k - q|) 
V 2 {Hq)F 2 (-k, k - q)P L (\k - q|) - ^ P L (q ) 


(A.12) 


In the last step we have absorbed the k —> 0 contribution of Cq\q E into the cr 2 s PL{k) term. 
Note that evaluation of the C^Q E (k) term in A.ll is more straightforward in the original 
coordinates of eq. A. 10. 
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A.1.3 The (5o5|<5) term 
Expanding the term to one-loop order: 

(<5o<S|(5) = (<jW o<yW|J (2) > + 2(^ 2) o<5W|<y (1) ). (A.13) 

As with P 0 \g(k), the smoothing filter factors out: 

Poo\s(k) = S(k)P 00 \ 0 (k) = S(k) [A 00 |o(fc) + 2B 00 \ 0 (k)] i (A.14) 

where we have: 

A 0 o|o(fc) = 2 J ^^F 2 (k- q , q)P L (g)P L (|k-q|), 

B 00]0 (k) = 2P L (k) j ^F 2 (-k, -q )P L (q) = 2 P L (k ). (A.15) 

Bringing everything together, we find: 

P 00l s(k) = S(k) ™a 2 P L (k) + 2 f ^F 2 (k-q, q )P L ( q )P L (\k-ql) . (A.16) 

A.1.4 The (S o 5|<5 o ( fiS)) c term 
Expanding the term to one-loop order: 

(5 o o (/+<*)) = (<$W o ^l^ 1 ) o (/+<5 (1) )). (A.17) 


This leads to the power spectrum: 

-Poo|oe(^) = 2 J j^M- q )P L (q)P L (\k- q \). 

Using eq. B.2 yields: 

P ooloE (k) = 2 S in 2 (0 k ) J ^V 2 (fX q )S(q)P L (q)(P L ( |k-q|)). 

Note that the k —> 0 limit is zero due to the quadrupolar dependence on fjL q . 


d 3 q 


A. 1.5 Combining terms 

Putting all the terms from eq. A.3 together yields: 


CO 

P g E{k,0 k ) =A(z) sm z (0 k ) x \bxPQ\sik) + b\—a 2 s P L {k) 

d 3 q 


+ 26? 


(27t) 3 


+ 2b(P L (k) J 

d 3 q 


V 2 {dq)F 2 {k - q, q )S(q)P L (q)P L {\k - q|) 
d 3 q 


iS(q) 


(A.18) 


(A.19) 


k,k-q)P L (|k-q|) - — -P L {q) 


+ 2b 2 


+ 26i6 2 


( 2 < 


(27r) 3 

3 F 2 (k - q, q )S(q)P L (q)P L (\k - q|) 


d 3 q 

(27r) ; 


V 2 ^ q )S(q)P L (q)(P L (\k- q \))+0(P 3 L ) 


(A.20) 
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where we have absorbed the divergent a l 2 contribution into the definition of b\ with the bias 
renormalization: 

68 

h -> h + —M 2 - (A.21) 

Combined with the earlier bias renormalization (b\ —> b\ + 3 & 30 -2 ), we have: 

bi —>► b \ + ^363 + —62^ v 2 , (A. 22 ) 

which is consistent with previous results (e.g. [47]). 


A.2 Calculating Pee and Per 

Using our freedom to choose k y = 0, the E- and .B-mode autocorrelation power spectra are 
given by: 


<7+ 1 7+) =A\z) x [(1 - 2b 2 a 2 )fl(5\5) + 2b 1 f + (5\6o(f + 5)) 

+ 2b 2 f + (8\8 o <5 o (/+$)> + b 2 (S o (f + 8)\8 o (/ + S)) + ML) 


<7 X |7x) =A 2 (z) x b\(5o (f x 6)\So (f x S)) + ML) 


(A.23) 


where we have used fx(k y = 0) = 0. Separating into disconnected and tightly coupled 
correlators, we find: 

( 7 +I 7 +} =A 2 (z) x [fl(5\5) + 2hU(S\S o (/+*)> + bl(6 O (f + 6)\6 o (/+$)> + 0(8f in ) 
<7x|7x) =^ 2 (*) x \bK 80 (f x 8 )\ 8 o (f x 8 )) + 0(SL) 


The (<5|<5) and (5\5o(f + 5)} terms have been described above. As expected, the leading B-mode 
contribution appears at one-loop and is due to weighting by the shape tracer density. 


A.2.1 The (5 o (fi8)\5 o ( faS )) term 

At one-loop, the only contribution to this term is (5^ o (fid^)\8^ o The B-mode 

correlation is: 

M:|o£( k ) = j + PL(q)PL(q2 ), (A.24) 

where = k — q. P 0 s| 0 B(k) is given by the equivalent expression with /e -» fs- Using 
eq. B.3, and making the Limber approximation (Qy. = vr/2): 


P ™\0E(k) -/^ 


3 - 14 /j 2 + 19 n\ 


) 


P L (q)S 2 (q) [P L {q 2 )-2P L {q)} 


8an q (3fi 2 - 1) + 4(—1 + 3 n 2 ) + a 2 (3 - 14 /i 2 + 19^) 


l 8(1 + a 2 — 2a n q ) 

P 0 B\0B(k) = J ^3 {2^(l-^))P L { q )S\q) (p L ( q2 )-2P L ( q )j 


PL(q)P L (q 2 )S(q)S(q 2 ) 


+ 




(A.25) 
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where a = q/k. Smoothing affects these terms on large-scales because of a zero-lag (scale- 
independent) contribution, which comes from a convolution over all scales (including the 
smoothed small scales). 


A.2.2 Combining terms 

Combining all contributions to the IA autocorrelation at one loop, we have: 

P EE ( k) = A\k) [sm 4 (6 k )S 2 (k)P 0l0 (k) + 2 b 1 sm\e k )S(k)P 0l0E (k ) 
+blP 0E \0E(k,e k ) + 0(P E )] , 

Pbb( k) = A 2 (k) [i b 2 P 0Bl0B (k,e k ) + 0(Pl)\ . (A.26) 


B Simplifying the angular integration of IA operators 


As seen in appendix A, the SPT expansion leads to convolutions of the linear power spectrum, 
similar in form to: 


p ( k ) ~ ^3 J d 3 q /*(q) F ( k - q>q)- p (9)^(l k - q|), (B- 1 ) 


where F denotes a gravity kernel, and fi is the relevant IA operator. We have the freedom 
to choose an arbitrary coordinate system for the wavevector q. If we rotate to a coordinate 
system in which k is the polar axis, the only azimuthal dependence in the integrand is in 
fi( q). We can then analytically evaluate the 4> q part of the integral. In this appendix we 
show the results of the azimuthal integrals of the relevant IA operators. We define q2 = k —q 
and use the freedom to set k y = 0 in the final results. 

For convolution integrals with only one IA operator, only /+ is relevant. 


1 

27T 


2tt 

J d(j) q f + { q) 
0 


1 

27T 


2tt 

J #<?/+( q2) 

0 


7,(3 k-q - 1 )f(k) = V 2 (n q )sm 2 (6 k ), 


1 + a 2 V 2 (n q ) - 2 a/iq 
1 -fa 2 — 2a n q 


sin 2 (6» fe ), 


(B.2) 


where Vi are the Legendre polynomials, a = q/k, fj, q = q ■ k = cos 6 q . 

For convolution integrals with two IA operators, the ^-dependence does not factor out 
of resulting expressions. For simplicity, we make the Limber approximation, in which only 
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transverse modes {Of. = 7r/2) contribute: 


2 7T 


^ f # 9 /+(q)/+(q) = ^(3 - 14/x 2 + 19/4), 


0 

27T 


# g /+(q)/+(q2) = 


1 

27r 

4 /#,/ + (q 2 )/ + (q 2 )= (1 + _ 2 )2 

0 


0 

2?r 


8a/i q (l - 3/x 2 ) + 4(3/4 - 1) + a 2 (3 - 14/4 + 19/4) 
8(1 + a 2 — 2 afi q ) 

1 — 4a/i g + 2a 3 /r g (l — 3/u 2 ) 


+ 


8a 2 (-l + 7/x 2 ) + a 4 (3 - 14/x 2 + 19/x 4 ) 
8(1 + a 2 — 2 ap. q ) 2 


2tt 


^ J # g /x(q)/x(q) = 2/r 2 (i-/x 2 ), 

0 

2tt 9 

1 , 2a/i,(/i,-l)(l-a/i 9 ) 

2iF J # ’ / x (q)/ '< (q2) = 1 + a 2 — 

0 


2?r 


[ #g/x(q2)/x(q2) = 


2a 2 (l -/x 2 )(l -a/Xg) s 
(1 + a 2 — 2 a/tq) 2 


(B.3) 


Note that because the IA operators are even, in the A; —>• 0 limit, these expressions are 
unchanged under the exchange q -f-t q 2 • 
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